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ABSTRACT 

Particle acceleration at astrophysical shocks may be very efficient if magnetic 
scattering is self-generated by the same particles. This nonlinear process adds 
to the nonlinear modification of the shock due to the dynamical reaction of 
the accelerated particles on the shock. Building on a previous general solution 
of the problem of par ticle acceleration with arbitrary diffusion coefficients 



lAmato fc Blasi 



(2005|)), we present here the first semi-analytical calculation 
of particle acceleration with both effects taken into account at the same time: 
charged particles are accelerated in the background of Alfven waves that they 
generate due to the streaming instability, and modify the dynamics of the 
plasma in the shock vicinity. 
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1 INTRODUCTION 



Soon a fter the pioneering papers by iKrvmskiil ()l977f l ; iBlandford fc Ostrikeii (|l978^ ; 



Bell 



|1978al lbfl. introducing the test particle theory of particle acceleration at coUisionless shocks, 
it became clear that the dynamical reaction of the accelerated particles on the plasmas in- 
volved in the shock formation may not be negligible. It is now clear that such reaction may 
in fact make shocks efficient accelerators and change quite drastically the predictions of the 
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test particle theory. The main consequences of the shock modification induced by the accel- 
erated particles can be summarized as follows: 1) a precursor, consisting in a gradual braking 
of the upstream fiuid, is created; 2) particles with different momenta feel different effective 
compression factors, which reflects in the fact that the spectrum of accelerated particles is 
no longer a power law, but rather a concave spectrum, as hard as p~'^-'^ at high momenta; 3) 
the shock becomes less efficient in heating the background plasma, so that the temperature 
of the downstream gas is expected to be lower than pr edicted throu g h the usual Rankine 
Hugon i ot relations at an unrn o dified shock frori t (see iDrurvl ()l983( ): iBlandford fc Eichle: 
( 19871 ): I.Tones fc EllisonI ()l991l ): iMalkov fc Drurvl ()20nil ) for reviews on different aspects of 
the subject). The reaction of the accelerated particles has been calculate d within differ- 



ent approaches: th e so-called two-fluid models f Prurv fc Vo 
els JMalkovl (|l997|); 



Malkov. Diamond &: Volk 



20001); 



Blasi 



jl980 



(2002 



appr o aches , both Monte Carlo and other s i mulat io n procedures 



198l|)), kinetic mod- 



2004^) and nu r nerica l 



Kang .Tones 



hson 


Vlobnis 


I997I 


2005 





.Tones Ellison 
fc PaschmannI (1990 ): lEUison. Baring fc .TonesI ()l995 



122a); 



Kang. .Tones fc Gieselerl ()2002| )). In most of these calculations. 



the diffusion properties of the plasma upstream and downstream are provided as an input to 
the problem. This also results in flxing the value of the maximum momentum of the accel- 
erated particles. However, one of the well known and most disturbing problems associated 
with the mechanism of particle acceleration at shock f ronts is that a substantial a mount 
of magnetic scattering of the particles is required (e.g. Lagage fc Cesarsky |l988al )bl)). In 
the absence of it, the maximum energy of the a ccelerated pa rt icles i s excee dingly low and 
uninteresting for astrophysical applications (e.g. iBlasil (|2005| )). iBelll (|l978al ) proposed that 
the streaming instability of cosmic rays could be responsible for the generation of pertur- 
bations in the magnetic fleld of an amplitude necessary to p rovide pitch angle sca.ttering 
(and therefore spatial diffusion) of the accelerated particles. iLagage fc Cesarskvl ()l988b( ) 
used this argument to estimate the maximum energy of particles accelerated at shocks in 
supernova remnants. In all previous works either the shock was considered unmodifled, or 
the diffusion coefficients were fix ed a priori, because a comprehensive theory of particle ac- 



celeration was missing. Recently, 



Amato &: Blasil (j2005|) have found a general exact solution 



of the system of equations describing the diffusion-convection of accelerated particles, and 
the dynamics and thermodynamics of plasmas in the shock region, for an arbitrary choice of 
the spatial and momentum dep endence of the d i ffusio n coefficient. In the present paper we 
use the formalism proposed by lAmato fc Blasil (|2005l ) and combine it with calculations of 
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the perturbations created through streaming instabihty, so that the diffusion coefficient, as 
a function of spatial location and momentum, is determined from the spectrum and spatial 
distribution of the accelerated particles. This provides the first combined description of the 
process of particle acceleration at collisionless shocks in the presence of particle reaction and 
wave generation. In this approach, the spectrum of accelerated particles, their distribution 

in the upstream plasma and the diffusion coefficient are outputs of the pro blem. 

Th e paper is organized as follows: in Sec. |21we summarize the findings of lAmato fc Blasi 



( 20051 ) . In Sec. El we illustrate our treatment of the streaming instability and determine a 
relation between the power spectrum of magnetic fluctuations and the diffusion coefficient 
upstream. In Sec. El we describe the results of our calculations. In Sec. El we shortly discuss 
how the results presented in the previous section change when the effects of turbulent heating 
are taken into account. We conclude in Sec. El 



2 CALCULATIONS OF THE SPECTRUM FOR ARBITRARY DIFFUSION 
COEFFICIENT 



In this section we briefly summarize the mathematical procedure proposed by 



Amato fc Blasi 



( 20051 ) to calculate the spectrum and spatial distribution of particles accelerated at astro- 



physical shocks, and their dynamical reaction on the shock stru cture, for an arbitrary dif- 



fusion coefficient D{x,p). The reader is referred to the paper bv lAmato fc Blasil ()2005r ) for 
more details. 

The equation for the conservation of the momentum between upstream infinity and a 
point X in the upstream region can be written as: 

where ^c{x) = Pcr{x)/ pou^ and U{x) = u{x)/uo and we used conservation of mass pqUq = 
p{x)u{x) (here po and Uq refer to the density and plasma velocity at upstream infinity, while 
p(x) and u{x) are the density and velocity at the location x upstream. Mq is the sonic Mach 
number at upstream infinity). 

The pressure in the form of accelerated particles is defined as 

\ PPmax 

Pcr{x) = - dp47cp'v{p)f{x,p), (2) 



3 



Vifij 



and f{x,p) is the distribution function of accelerated particles. Here Pinj and Pmax are the 
injection and maximum momentum. The function / vanishes at upstream infinity, which 
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implies that there are no cosmic rays infinitely distant from the shock in the upstream 
region ^ . The distribution function satisfies the following transport equation in the reference 
frame of the shock: 



d_ 

dx 



d 

D{x,p)—f{x,p) 



dx 3 \dx , 



(3) 



The X axis is oriented from upstream infinity (x = — oo) to downstream infinity {x = 
+oo), with the shock located at a; = 0. The injection is introduced here through the function 
Q{x,p). The diffusion properties are described by the arbitrary function D{x,p), depending 



on 



3oth momentum and s pace ^. 



Amato &: Blasil (|2005| ) showed that an excellent approximation to the solution f{x,p) 
has the form 



f{x,p) = fo{p) exp 



(lip) f° u{x') 



3 Jx D{x',p) 

where fo{p) = f{x = 0,p) is the cosmic rays' distribution function at the shock and q{p) 
_ d.infoip) -^g local slope in momentum space 



(4) 



dlnp 

The function fo{p) can be written in a very general way as found by iBlasil (|2002| ): 

SRtot \ mo f r ^P' ^RtotUp{p' 



\RtotUp{p) - 1 j 47rp3^^. """"P I p' RtotUpip') - 1 

Here we introduced the function Up{p) = Up/uo, with 
1 ro 



u 



Ui - 



dx{du/dx)f{x,p) 



(5) 



(6) 



where Ui is the fluid velocity immediately upstream (at x = 0^). We used Q{x,p) = 
r?ngaa,itti ^^^ -Pmj)5(x), with Ugas,! = uoRtot/Rsub the gas density immediately upstream 
{x = 0~) and t] the fraction of the particles crossing the shock which are going to take part 
in the acceleration process. In the expressions above we also introduced the two quantities 
Rsub = Ui/u2 (compression factor at the subshock) and Rtot = uq/u2 (total compression fac- 



^ This assumption implies that we are not considering any reacceleration of pre-existing seed particles. 

^ Since we will be using this equation in Sect.|3]for the case in which diffusion is due to a strongly amplified turbulent magnetic 
field, a few comments are in order: rigorously, this equation describes the isotropic part of the distribution function, and as 
long as the quasi-linear theory holds, the anisotropic part is expected to represent a small perturbation. It is not clear how the 
equation would generalize to the strongly non-linear case, though it may be reasonable to assume that the anisotropy remains 
rather small as long as the Alfven speed in the perturbed field is negligible compared with the fluid speed. 

^ In writing Eq. EJin this form we are neglecting the velocity of the scattering centers with respect to the fluid velocity 
upstream. This is always a good approximation for the cases considered in this paper. 
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tor). If the heating of the upstream plasma takes place only due to adiab atic compre ssion, 
the two compression factors are related through the following expression (jBlasil ()2n02l )): 



tot 



(7. + ml, - (7. - 1)^; 



79 + 11 

sub 



(7) 



where Mq is the Mach number of the fluid at upstream infinity and 7g is the ratio of specific 
heats for the fluid. The parameter rj in Eq. El contains the very important information about 
the iniection of particles from the thermal bath. We adopt here the recipe proposed by 
Blasi. Gabici fc Vannonil (|2005i ) that allows us to relate rj to the compression factor at the 



subshock as: 
4 



7] 



37rV2 



{Rsuh - l)re" 



(8) 



Here ^ is a parameter that identifies the injection momentum as a multiple of the momentum 

of the thermal particles in the downstream section {pmj = C,Pth,2)- The latter is an output of 

the non linear calculation, since we solve exactly the modified Rankine-Hugoniot relations 

together with the cosmic rays' transport equation. For the numerical calculations that follow 

we always use ^ = 3.5, that corresponds to a fraction of order 10~^ of the particles crossing 

the shock to be injected in the accelerator. 

In terms of the distribution function (Eq. 0}, we can also write the normalized pressure 

in accelerated particles as: 
Air 



dp p v{p)fo{p) exp 



Xp{x',p) 



3D(p,x) 
q{p)uo 



3po«o 

where for simplicity we introduced Xp{x,p) 

By differentiating Eq. IHl with respect to x we obtain 

dx 

where 



(9) 



X{x)^c{x)U{x), 



A(x) =< 1/xp >^ 





\ fOj^l U{x') 1 
[ Jx "-^ Xp(x',p)\ 


IpZr P^^iP)fo{p) exp 


•'X Xp(x,p) 





(10) 



(11) 



and U (x) is expressed as a function of C,c{x) through Eq. [TJ 

Finally, after integration by parts of Eq. IHl one is able to express Up{p) in terms of an 
integration involving U (x) alone: 

1 



Upip) 



dx U{x) 



Xp(x, p) 



exp 



'dx' 



(12) 



Xp{x',p)_ 

which allows one to easily calculate fo{p) through Eq. El 

Eqs. H and ^1 can be solved by iteration in the following way: for a fixed value of the 
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compression factor at the subshock, Rsub-, the value of the dimensionless velocity at the shock 
is calculated as ?7(0) = Rsub/Rtot- The corresponding pressure in the form of accelerated 
particles is given by Eq. [T] as ^c(O) = 1 + - - ("^) '^'^^^ ^^^^ ^ 
boundary condition for Eq. ^1 where the functions U{x) and X{x) (and therefore fo{p)) on 
the right hand side at the fc*^ step of iteration are taken as the functions at the step {k — 1). 
In this way the solution of Eq. ^ at the step k is simply 



(a:)=ec(0)exp 







(13) 



with the correct limits when x ^ and x — * — oo. At each step of iteration the functions 
U{x), fo{p), X{x) are recalculated (through Eq. [H Eqs. IT^ and 131 and Eq. [TTl respectively), 
until convergence is reached. The solution of this set of equations, however, is also a solution 
of our physical problem only if the pressure in the form of accelerated particles as given by 
Eq. n coincides with that calculated by using the final fo{p) in Eq. IHl This occurs only for 
one specific value of Rsub, which fully determines the solution of our problem for an arbitrary 
diffusion coefficient as a function of location and momentum. 



3 SELF-GENERATED TURBULENCE AND PARTICLE DIFFUSION 



The streaming of cosmic rays at super- Alfveni c speed indu ces a streaming instability, which 
has been discussed in previous literature fe.g. iBelll (1978a)). 

Let us define J-'{x, k) as the energy density per unit logarithmic band width of waves 
with wave-number k. Neglecting the damping, and assuming a steady state, the following 
relation holds (see e.g. Lagage & Cesarsky 1983): 
(9jF(x, k) 



u- 



dx 



a{x, k)T{x, k) 



(14) 



4 df{x,p) 
p V 



dx 



(15) 



where u = u{x) is the fiuid velocity upstream of the shock and a is the growth rate of waves 
with given wavenumber k, which can be related to the distribution function of the resonant 
cosmic rays, f{x,p{k)), through: 
4 vr Va 

In Eq. El V and p are the particle velocity and momentum respectively, and the latter is 
related to the wave number k through the resonance condition p{k) = eB/kmc, Um is the 
energy density of the background magnetic field Bq {Um = Bq/Stt), while va is the local 
Alfven velocity. 



p=p{k) 
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All these expressions have been actually obtained for shocks that are not modified by 
the dynamical reaction of cosmic rays. In principle, the Fourier analysis used to obtain 
the previous expressions and in fact used to reach the conclusion that there are unstable 
modes, is not formally applicable, since all these calculations assume that the background 
quantities (the fiuid velocity u in particular) are spatially constant. However, provided that 
1/k remains much smaller than the spatial extension of the precursor, the conclusions are, 
in first approximation, still applicable. Clearly this condition is broken by definition at the 
maximum momentum Pmax at least in those cases in which this is determined by the finite 
size of the accelerator rather than by energy losses. Special care should be taken of the fact 
that all quantities involved in the equations above depend on the location in the presursor. 

It follows that for a cosmic ray modified shock, va is not spatially constant since both the 
upstream plasma density, p, and, in general, the background magnetic field, i?o, are space 
dependent. However all previous calculations apply to the case of a parallel shock, for which 
the strength of the background magnetic field Bq can be taken to be constant, since there 
is no adiabatic compression of the magnetic field lines. 

Using the equation for conservation of mass p{x) = po/U{x), we can therefore write the 
local Alfven velocity as: 

^^^> - Tfe"'^'"^ • ^''^ 

An additional warning should be issued in that Eg. IT^ neglects the adiabatic compression 
of waves in the shock precursor: this refiects in the absence of terms proportional to the 
gradient of the velocity field. Unfortunately, to our knowledge, discussions of this problem 
in the literature are limited to integrated quantities (e.g. total energy density and pressure 
of the waves) while a description of the behaviour of the modes with different wave-numbers 
is more complex. In fact, in principle even the concept of modes with given k becomes ill 
defined in a background which has spatial gradients of the quantities to b e pertu rbed. 



Once J-'{x, k) is known, the diffusion coefficient is known in turn (jBelll (|l978a|)): 



D{x,v) = -'^. (17) 

From the latter equation, where stands for the Larmor radius of particles of momentum 
p, it is clear that the diffusion coefficient tends to Bohm's expression for JF — > 1. On the 
other hand, it is also straightforward to check what the expected saturation level for the 
overall energy density of the perturbed magnetic field is. If we define 
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6B^ _ B^ rdk 



87r 87r y A; ^ ^ ' ^ ^ 

from Eq. IT^and Eq. ^Jwe see that: 
dl Bl rdk 

""Tx = 8^yT"-^(^'^)= 

= j dp v{p) f{x,p) = (19) 

dx 

Integration of the latter equation is straightforward when non-hnear effects on the fluid are 
neglected so that u and va are both spatially constant. One obtains 6B'^/8tt = {va/u)Pcr, 
or, in terms of amplification of the ambient magnetic field: 

(fX = 2MA^, (20) 

with Ma = Uq/va the Alfvenic Mach number. 

It is worth stressing that for Pcr/ Po^l ~ 1 and Ma ^ 1, the predicted amplification 
of the magnetic field exceeds unity. In fact, this result was initially obtained in the context 
of the so-called quasi-linear theory, therefore it should be taken with caution and checked 
versus numerical calculations of the non-linear phase of amplification of the waves. It seems 
clear, however, that the growth may well enter this non-linear regime and lead to turbulent 
fields in the shock vicinity that exceed the pre-existing background magnetic field. 

Let us now go back to Eqs. [T^HTT] with the aim of recasting the relation between the 
diffusion coefficient and the cosmic ray distribution function in a more compact form. Using 
Eqs. [T^HTB] . we can write: 

^{x,P) = — , (21) 

3 PO Uo f AO 

where vao = Bq / y/4:7ipo is the Alfven speed at upstream infinity and 

With this definition of $, from Eq. El we obtain: 

p.. X 3 no VAO Uo 

Di-^P) = ^2DBo^^^^—-, (23) 

with Dbo = nipC^ /3eBo a constant. 

It is important to notice that since the constant Dbo is inversely proportional to the 
strength of the background field Bo, and the Alfven speed vao is proportional to Bo, the 
diffusion coefficient in Eq. 1221 turns out to be independent of Bo- This result holds only 
within the context of quasi-linear theory. Even in the context of a quasi-linear theory of the 
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development of magnetic perturbations, a dependence on Bq could be introduced through 
the quantity $, which is affected by the laws of conservation of momentum and energy in 
the precursor. However, in the cases of interest for us we will show below that these effects 
are fully negligible. 

4 SPECTRA OF THE ACCELERATED PARTICLES AND 
SELF-GENERATED DIFFUSION COEFFICIENT 

Following the mathematical procedure outlined in Sec. |21 and in Sec. 01 we are able to de- 
termine self-consistently the spectrum of accelerated particles and the diffusion coefficient. 
Within the obvious limitation of using quasi-linear theory to calculate the diffusion coefficient 
for the non-linear case, this is the first attempt at determining the space and momentum 
dependence of the diffusion coefficient together with the spectrum of accelerated particles. 
While in a time-dependent approach to the problem it would be possible to estimate the 
maximum energy in a self-consistent way, here we assume for simplicity that the maximum 
momentum is a given parameter. We chose to carry out the calculations presented in the 
following for Pmax = lO^mc. 

The spectra of the accelerated particles for Mach numbers at upstream infinity ranging 
from Mq = 4 to Mq = 200 are shown in Fig. ^ for a background magnetic field at upstream 
infinity Bq = IfiG. As stressed above the result is however expected and actually found to 
be independent of the strength of the background magnetic field. In the bottom part of the 
same figure we plot the slope of the spectrum as a function of momentum. 

It is evident that for low Mach numbers and at given Pmax the modification of the shock 
due to the reaction of the accelerated particles is small (see for instance the case Mq = 4). 
For the strongly modified case (e.g. Mq = 200) the asymptotic spectrum of the accelerated 
particles is very flat, tending to p~°' with a = 3.1 — 3.2 for p — > Pmax- The momentum at 
which the spectrum becomes flatter than the prediction of linear theory, depends on the 
level of shock modification: it is higher (10—20 mc ) for relatively low Mach numbers (namely 
weaker modification) and approaches a few GeV for high Mach numbers and large shock 
modification. The asymptotic spectrum is reached at p/mc > 10^. These effects might be 
important in the perspective of reconciling the concave shape of the instantaneous spectra of 
accelerated particles with observations of the diffuse spectrum of cosmic rays in the Galaxy. 
Most measurements, mainly related to the abundance of light elements are in fact limited to 
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p / mc 



Figure 1. Spectrum and slope at the shock location as functions of energy for pmax = 10^ mc and magnetic field at upstream 
infinity Bo = l/^G. The curves refer to Mach numbers at upstream infinity ranging from Mq = 4 to Mq = 200: dotted for 
Mo = 4, dashed for Mq = 10, dot-dashed for Mq = 50, solid for Mo = 100 and dot-dot-dashed for Mo = 200. 

relatively low energies, where the spectra predicted in this paper are compatible with power 
laws softer than p~^. Serious work aimed at predicting the actual spectrum of cosmic rays 
escaping the sources is urgently needed but still missing in the context of non-linear theories 
of particle acceleration at shocks. 

The diffusion coefficient associated with the self-generated waves is given by Eq. |221 We 
plot this diffusion coefficient at the shock location in Fig. |21 for Mach numbers Mq = 10 
(dashed lines) and Mq = 100 (solid lines). We fix Bq = IfiG, but as stressed in the pre- 
vious section, the diffusion coefficient obtained within the quasi-linear theory of magnetic 
perturbations is independent of Bq. For comparison, we also plot the corresponding Bohm 
diffusion coefficient Dsip) oc v{p)p in the unperturbed magnetic field Bq, for Bq = l/iG 
and Bq = lO^G. The comparison strikingly shows that for most momenta of the accelerated 
particles the diffusion takes place at super-Bohm rates (namely the diffusion is slower than 
predicted by the Bohm coefficient in the unperturbed magnetic field, as could be expected). 
Moreover, the difference between the self-generated diffusion coefficient and the Bohm coeffi- 
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cient increases at the highest momenta, which might suggest that somewhat higher energies 
could be achieved if the self-generated turbulence were taken into account. In this respect 
it is also important to notice that, as one might expect, the more modified the shock, the 
slower the diffusion. 

It is worth keeping in mind that the diffusion coefficient in the amplified magnetic 
field, as obtained through our calculations, remains larger than the Bohm value in the 



same field. The latter is in 



act considered as a sort of lower limit to the diffusion rate 



( Casse. Lemoine &: Pelletied (|2002i )) even in the case of strong turbulence. The only region 
in momentum space where this condition may be violated in our calculations is very close 
to the maximum momentum Pmax- If is clear however that a realistic determination of the 
diffusion coefficient cannot be achieved in the conte xt of quasi-linear theorv an d that even 



numerical approaches to diffusion, such as those of ICasse. Lemoine fc Pelletieii ()2nn2f ) can 
only suggest a general trend as long as the turbulent structure of the magnetic field is pre- 
defined rather than determined by the diffusing particles themselves. In this sense, the limit 
at the Bohm value in the amplified field should also be taken with caution. 

As stressed above, the fact that the diffusion coefficient is smaller than the Bohm co- 
efficient in the background field is the consequence of the fact that the fluctuations in the 
magnetic field become strongly non linear, namely 6B'^/Bq ^ 1, at least close to the shock 
surface. In fact we find that 6B/Bq at x = is exactly as predicted by Eq. 1201 In these 
conditions it is important to check that the dynamical role of the turbulent magnetic field 
remains small. In Fig. IHlwe plot SB'^/Stt normalized to PqUq (top panels) and the cosmic 
ray normalized pressure C,c{x) and velocity U{x) (bottom panels). The curves refer to Mach 
number Mq = 10 (dashed lines) and Mq = 100 (solid lines). The plots on the left (right) are 
obtained for Bq = lO^G {Bq = IjiG). The x-coordinate is in units of = — -DB(pmax)/wo, 
where Db{p) stands for the Bohm diffusion coefficient appropriate to the considered value 

of Bq. 

The highest values of SB"^ /SttpqUq, reached close to the shock front, are of the order 
of 10^^ — 10^^, confirming that even in the extreme non linear cases the dynamical effect 
of the magnetic field remains unimporant. This result serves as a justification a posteriori 
that we could neglect the pressure of the waves and their energy flux in the equations 
of conservation of momentum and energy respectively. This result is very specific of the 
resonant channel of production of Alfven waves, and is very likely not correct in the case of 
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u 



Q_ 
Q 




10- 



10° 



10^ 



10^ 



(p/mc) 



Figure 2. The self-generated diffusion coefficient at the shock location a; = 0~ as a function of the particle momentum for 
Mach numbers Mo = 10 (dotted line), Mo = 100 (dashed line) and Mo = 200 (solid line). Also plotted is the Bohm diffusion 
coefficient corresponding to Bo = l^^G (solid line with triangles) and Bo = 10/iG (solid line with diamonds). The y-axis is in 
units of cm^s~^. 




=1/iG 




10-10 10-8 10-6 10-* 10-2 10° 10-'0 10-8 10-6 IQ-" 10-2 10° 

x/x. x/x. 

Figure 3. Top panels: the energy density in magnetic field fluctuations 5B^ /Sw normalized to the fluid ram pressure pqUq 
at upstream infinity. Bottom panels: the eosime ray pressure normalized to PoMq, 5c (thick curves), is plotted together with 
the normalized velocity U (thin curves). All functions are plotted versus spatial location, vifith the x-eoordinate in units of 
I* = —Dg(pms.yi)/uQ, where Db(p) stands for the appropriate Bohm diffusion coefficient. The left and right panels refer to 
different strengths of the background magnetic field Bu, as specified in each panel, while the different line-types correspond 
to different Mach numbers: dashed for Mo = 10 and solid for Mo = 100. In the upper panels we also plot for comparison a 
dot-dashed curve corresponding to SB = lOBo. 
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10-2 IQO 102 104 

p / mc 

Figure 4. Effects of using the self-generated diffusion coefficient on the particles' distribution functions. In the top and bottom 
panels we plot the particles' spectrum and spectrum slope respectively. The continuous curves are obtained for self-generated 
D{x,p): dashed curves are for Mg = 10 and solid for Mg = 100. The symbols represent the results obtained adopting the Bohm 
diffusion coefficient: diamonds are for Mq = 10 and filled circles are for Mq = 100. 



non- resonant scenarios, such as the one proposed bv iBelll (200J), where larger amphfications 
of the magnetic field could be achieved. 

The shape of the spectra of accelerated particles is affected in a sizeable way by the 
adoption of the self-generated diffusion coefficient: Fig. 0] illustrates this point. The continu- 
ous lines are for self-generated diffusion (dashed for Mq = 10 and solid for Mq = 100) while 
the symbols are for Bohm diffusion (diamonds for Mq = 10, almost perfectly superposed 
on the dashed curve, and filled circles for Mq = 100). While in the weakly modified cases 
the spectrum is basically independent of the assumed diffusion coefficient, in the fully non 
linear solution self-generated diffusion leads to steeper spectra at low momenta and harder 
spectra at high momenta. 



14 E. Amato and P. Blasi 

5 THE ROLE OF TURBULENT HEATING 

A major uncertainty in all types of calculations of the non linear particle acceleration at 
shock fronts is the effect of turbulent heating. This generic expression is used to refer to any 
process that may determine non-adiabat ic gas heating. The two best known examples of 
this type of pro c esses are Alfven heating (jMcKenzie fc Volkl ()l982l )) and acoustic instability 
( Drurv fc Falld (|l986| )). Both effects are however very hard to implement in a quantitative 
calculation: in the case of Alfven heating, the mechanism was originally introduced as a way 
to avoid the turbulent magnetic field to grow to non linear levels, while it is usually used 
even in those cases in which 6B/Bq ^ 1. 

Acoustic instability develops in the pressure gradient induced by cosmic rays in the 
prec ursor and resul t s in th e development of a train of shock waves that heat the background 
gas (jPrurv fc Falld (|l986r i). The analysis of the instability is carried out in the linear regime, 
therefore it is not easy to describe quantitatively the heating effect. 

In both cases the net effect is the non-adiabatic heating of the gas in the precursor, 
which results in the weakening of the precursor itself and in the reduction of the acceleration 
efficiency compared with the case in which the turbulent heating is not taken into account. 



of 



In order to illustra t e this effect, we adopt a phenomenological approach, similar to that 



Berezhko &: EUisonl ()l999| ). We stress that this approach, developed for the case of weak 



turbulence, is inadequate in principle for the case of interest here, where the turbulence 
can, in principle, become strong. We adopt it here only for illustration of the main physical 
effects. 

The approach consists in redefining the equat ion of state of the gas ta king into account 
the heating induced by the accelerated particles (|McKenzie &: Volkl (|l982l )): 



d 



vh{x) dPcR 



p{x) 



-79 



(24) 



dx ^ ' u{x) dx 

where u/vh = Mh is the local Mach number of the turbulence relevant for the heating (for 
instance vh = va for Alfvenic heating). After defining r(x) = Pg{x) / poul, we replace Eq. [T] 

in the set of equations presented in Sec. El with the two following equations: 

1 



Ux) + Uix)+T{x) 

U{x)-^^ 



1 + 



t{x) 



1 + 7,(79-1) 



M2 \1-U{x)^^- 



HO 



+ lJx) 



(25) 



(26) 



where the latter is obtained by rewriting Eq. |2^in terms of the normalized pressures and 
integrating between upstream infinity and a generic location x in the upstream medium. 
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after expressing ^d^) in terms of r(x) and U{x) through Eq. |2H1 We have assumed a spatial 
dependence of the turbulence characteristic velocity vh in the form of vh{x) = vho U{xY 
(with s = 1/2 in the case of Alfven heating, from Eg. Il(j|l. and used as a boundary condition 
t{—oo) = l/{pfgMQ). The term Xt-(x) appearing in Eq. 123 finally, is defined as: 

f/(a;)7.-i+^ 7^ dx' . (27) 

-oc CiX 

The only other changes induced by the inclusion of turbulent heating in our initial set of 
equations concern the relation between the compression ratios Rtot and Rsub (Eq. Ej) and the 
temperature jump between downstream and upstream infinity, that refiects on the minimum 
cosmic ray momentum pinj. In both cases the changes can be summarized in the appearence 
of a factor (1 + Ffj) with 



^ I Rsub \ '^^ 



Mho [lg + s [ V Rtot ) 
With this definition of Fu we find: 



+ X.(0) . (28) 



Rtot = M^'- 
and 



^sub 



2(1 + Fh) 



RtotV^~\, , ^ Aig + i)-{ig-m:ub 



(29) 



\Rsub/ {'^g + l) -{'^g-l)Rsub 

The usual results (adiabatic heating) are recovered when Mho/Mq —>■ oo. In spite of the 
apparent simplicity of these revised relations there are two complications arising in the 
solution of the system of equations. A minor difficulty is that the equation relating Rsub 
and Rtot now cannot be solved analytically due to the presence of the ratio Rsub/ Rtot in the 
definition of Fh- A more serious complication, instead, has to do with 2r, which requires 
the knowledge of the complete solution of the problem. However not even this is too severe 
a problem to overcome within the framework of an iterative method, although sometimes it 
results in an appreciable slowing down of the calculation. In fact, it can be seen a posteriori 
that 2r is always negligible compared to (1 — U'^'-'^^) / ijg + s). 

In order to show how our results for particle acceleration may be affected by the inclusion 
of turbulent heating we carried out the calculations for the case of Alfven heating, namely 
considering the Alfven velocity as the characteristic turbulence velocity, vh = v^, which 
also implies s = 1/2 in the equations above, according to Eq. ^1 Results obtained with and 
without inclusion of Alfven heating are shown in Fig. where for the background magnetic 
field we have assumed the largest of the two values so far considered, Bq = lO^G, with the 
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10-'0 10-'' 10-5 10-^ 10-2 IQO p/^^ 

x/x. 

Figure 5. The plots on the left show: in the upper panel, the ratio between the turbulent and background magnetic field as a 
function of space for two different values of the Mach number (10 and 100), with and without inclusion of the turbulent heating; 
in the lower panel the corresponding normalized cosmic ray pressure. The plots on the right show the particles' spectrum and 
slope in the same cases. The continuous curves correspond to cases when the turbulent heating is not taken into account: 
dashed for Mq = 10 and solid for Mq = 100. The symbols correspond to cases including the turbulen heating: diamonds for 
Mo = 10 and filled circles for Mq = 100. 



aim of minimizing the turbulence Mach number and hence maximizing the effects of the 
heating. 

In the left panel of Fig. El we show how the turbulent magnetic field strength and cosmic 
ray pressure are now reduced, for both a mildly (Mq = 10 diamonds versus dashed line) and 
a strongly (Mq = 100, filled circles versus solid line) modified case. In the strongly modified 
case we find that, in the vicinity of the shock, the turbulent magnetic field strength is 
decreased by 10%, while the cosmic ray pressure is decreased by 20%. Changes in both 
quantities are of order few % in the weakly modified case. In the right panel of the same 
figure we show, using the same notation for the different curves, how the particles' spectra 
are affected: while changes are negligible in the weakly modified case, for Mq = 100 the 
concavity of the spectrum is appreciably reduced, namely, the spectrum becomes harder 
toward the low energy end and softer at high energies. 

It is worth stressing once more that this way of including turbulent heating, that is used 
in many currect approaches to particle acceleration in supernova remnants, is far from self- 
consistent and the results should only be considered as an indication of a trend. Sometimes, 
in order to attempt a slightly more realistic approach, one substitutes the Alfven speed 
in the background field Bq with the corresponding quantity in the amplified field Bq + 5B. 
Needless to say that such an attempt, though justified by the complete lack of any non-linear 
theory of turbulent heating, is far from being realistic. 
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6 CONCLUSIONS 

We described the mathematical theory of particle acceleration at non-relativistic shock fronts 
with dynamical reaction of accelerated particles and self-generated scattering waves. The 
diffusion coefficient itself is an output of the calculations, though within the limitations 
imposed by the usage of quasi-linear theory applied to the case of potentially strong magnetic 
field amplification. The scattering in the upstream plasma is generated through streaming 
instability, as discussed extensively in previous literature. 

We determined the spectra of accelerated particles, their spatial distribution and the 
space dependence of the fluid velocity, pressure and temperature. The diffusion coefficient 
and the strength of the self-generated magnetic perturbations are also calculated, as a func- 
tion of the distance from the shock front in the precursor. We confirm the general finding 
that the spectra of accelerated particles are concave, an effect which is particularly evident 
for strongly modified shocks, namely for large Mach numbers of the moving fluid. However, 
the shape of the concavity is somewhat affected by the self-determined diffusion coefficient, 
as visible in Fig. |3] 

Having in mind the comparison between the predicted spectra at the sources and the 
observed cosmic ray spectrum at the Earth, it is worth reminding the reader that what can 
actually be measured is the combination of the diffusion time, the gas density along the 
trajectory (responsible for the spallation) and the injection spectrum. In order to infer some 
conclusions about the spectrum at the source, one has to make assumptions on the diffusion 
coefficient in the interstellar medium. In alternative it would be a precious step forward if 
we could measure unambiguously the spectrum of gamma rays generated by vr" decays close 
to the source itself, an evidence that unfortunately is still missing. 

The asymptotic slope of the spectra for p Pmax may be as flat as ~ 3.1 — 3.2, but 
this conclusion is not strongly affected by the fact that the diffusion coefficient is calculated 
self-consistently. 

The most striking new result of our calculations is the energy dependence of the diffusion 
coefficient and the strength of the amplified turbulent magnetic field. As could be expected, 
the diffusion coefficient is not Bohm-like, and the turbulent component of the magnetic 
field is amplified so efficiently that the diffusion coefficient is much smaller than the Bohm 
coefficient in the background magnetic field. This is especially true at the highest momenta, 
which leads to think that a full non-linear theory might predict higher values of the maximum 
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momentum than expected on naive grounds. Unfortunately a full, self-consistent calculation 
of Pmax for a strongly modified shock has never been carried out, the main difficulty being 
in accounting for the spatial dependence of all the quantities involved. 

When compared with the Bohm diffusion coefficient as calculated in the amplified mag- 
netic field, our diffusion coefficient remains always larger, with the possible exception of a 
narrow momentum region close to Pmax- 

While the calculation presented here is fully self-consistent in the determination of the 
shock modification due to the reaction of the accelerated particles, the part related to the 
amplification of the background field suffers from all the limitations related to the usage 
of quasi-linear theory for the streaming instability. This approach, initially developed for 
weakly amplified magnetic fields, is widely applied in the literature to situations that violate 
this condition. Unfortunately at the present time this is the only way we have to achieve a 
(at least partially) self-consistent picture of the process of particle acceleration at cosmic ray 
modified shocks with self-generated turbulence. This problem is in fact even more serious for 
those approaches that pred ict levels of niagneti c field ampli fication which are much higher 
than those found here (e.g. Bell fc Lucekl poniL Belli JsOOJ)). 

The high acceleration efficiencies obtained in the context of all approaches to particle 
acceleration at shocks are known to be reduced by the effect of turbulent heating. Any 
non-adiabatic heating of the gas in the precursor leads to reducing the energy channelled 
into non-thermal particles at the shock. This is a serious problem, because the effect of 
turbulent heating depends dramatically on the type of mechanism that is responsible for the 
heating: Alfven heating, often used in the literature, is only one of these mechanisms, and 
not necessarily the most efficient. For instance, the instability induced by the propagation 
of acoustic waves in the precursor was shown to lead to t he formation o f weak shocks in the 
precursor, which in turn heat the upstream plasma (e.g. 

These non-linear effects can hardly be taken into account in a credible way. Most notably, 
the phenomenological expressions proposed in the literature and used also in the present 
paper, have originally been proposed as mechanisms to reduce the amount of magnetic field 
amplification and remain in the context of small perturbations of the background magnetic 
field. However, as shown in Fig. 5, even with the Alfven heating taken into account, the 
magnetic field can be amplified by a factor in excess of ~ 10 with respect to the background 
field. This means that a fully non-linear theory of the turbulent heating is required in order 
to make fully reliable predictions. 



Drurvfc Fallel (|l98fi[l). 
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From the phenomenological point of view, the best evidences for both magnetic field am- 



phfication and effi cient particle acce 



fsee the reviews of 



Hi11a.' 



(200^ and 



eration com e from observations of supernova remnants 



Blasil ()2nn5r ) and references therein). In fact, it has been 



argued that the amount of field amplification required to expl ain the thickness of the X-ray 



bright rims in several remnants is of the order of ~ 200— 300/xG (|Volk.Berezhko. &: Ksenofontov 
( 20051 )). An important role in explaining thi s level of amp 



ferent versions of the streaming instability (|Bell &: Lucek 



ification cou 



(200l|); 



d be p l ayed by dif- 



Belll (j2004 . 



20051)), not 



requiring resonant interactions of particles and waves. A full non-linear theory including 
these effects will be described elsewhere (Amato and Blasi, in preparation). 
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